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We develop a collisional formalism adapted for the dynamics of ultracold dipolar particles in a 
confined geometry and in fields tilted relative to the confinement axis. Using tesseral harmonics 
instead of the usual spherical harmonics to expand the scattering wavefunction, we recover a good 
quantum number £ = ±1 which is conserved during the collision. We derive the general expression 
of the dipole-dipole interaction in this convenient basis set as a function of the polar and azimuthal 
angles of the fields. We apply the formalism to the collision of fermionic and bosonic polar KRb 
molecules in a tilted electric field and in a one-dimensional optical lattice. The presence of a 
tilted field drastically changes the magnitude of the reactive and inelastic rates as well as the 
inelastic threshold properties at vanishing collision energies. Setting an appropriate strength of the 
confinement for the fermionic system, we show that the ultracold particles can even further reduce 
their kinetic energy by inelastic excitation to higher states of the confinement trap. 


I. INTRODUCTION 

The field of ultracold gases composed of dipolar par¬ 
ticles hasgenerated tremendous interest during the past 
years |l|-|6|] . One major goal is to shape at will the quan¬ 
tum properties of ultracold gases using the high degree of 
controllability available in experiments. Different kinds 
of dipolar particles are concerned by the quest for man¬ 
ifestation of dipole-induced effects. The first category 
of interest contains particles with electric dipoles such 
as ground state molecules of KRb @, H| , RbCs [§, |Toj] , 
NaK [Tlj and many others experimentally under way... 
The second category includes particles with magnetic 
dipoles such as Cr fl2Ul4{ . Dy [lEIIl], Er [13,0] atoms 
and Er 2 molecules [Kj... The third category consists of 
particles with both electric and magnetic dipoles such as 
molecules of OH 0, SrF [H| , YO [H, E3 , RbSr [H]... 
All these particles can be manipulated by different con¬ 
figurations of electric and/or magnetic fields. They can 
also be loaded in optical lattices of different dimensions 
such as one dimensional (ID) lattices [25|-|29j], two dimen¬ 
sional (2D) lattices [30, [3l[ and three dimensional (3D) 
lattices |32|. Due to the wide and numerous domains of 
application of dipolar particles UK , it is therefore impor¬ 
tant to understand how to control the interactions and 
the collisional properties of these particles under such 
configurations of fields and lattices. For example, it has 
been shown that chemical reactivity of molecules can be 
suppressed by an electric field in a confined ID optical 
lattice m or by selecting a particular electric field and 
appropriate quantum states of the molecules in a non- 
confined space (34j . 

In this study, we investigate two-body collisions in a 
confining ID optical lattice in electric and/or magnetic 
fields tilted with respect to the £ and x axis, as illus¬ 
trated in Fig. [T] We consider here “classical” dipoles 
aligned along the field, for which the angular internal 
structure of the particles (rotational angular momentum 
for the electric dipoles and electronic angular momentum 
for the magnetic dipoles) is not taken into account. We 
choose an effective value d or /r of the electric or mag¬ 



FIG. 1: (Color online) Electric and magnetic dipolar collisions 
in a ID confinement (pancakes-shaped lattices) and arbitrary 
tilted fields relative to the a and x axis. We consider “classi¬ 
cal” dipoles aligned with the fields (see text). 


netic dipole moment, corresponding to their expectation 
value along the direction of the electric or magnetic field. 
For particles without angular internal structure, the total 
angular momentum J of the two-body colliding system 
reduces to the orbital angular momentum L associated 
with the quantum number l , which is conserved in free 
space. When a field is applied parallel to the quantiza¬ 
tion axis z, it is not conserved anymore but its projec¬ 
tion L z associated with the quantum number mi still is. 
When the field is tilted and no more parallel with respect 
to the quantization axis, its projection is not conserved 
anymore. The scattering problem becomes challenging as 
l, mi are all mixed. We show in this study that it is still 
possible to define a good quantum number, provided that 
the collision is described using the so-called “tesseral har¬ 
monics” [35I instead of the standard spherical harmonics 
for the partial wave expansion of the scattering wavefunc¬ 
tion. The problem is thus split in two sub-problems of 
smaller size. As an example of a dipolar system, we study 
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fermionic and bosonic KRb + KRb collisions in a tilted 
electric field. The present theoretical formalism was also 
successfully used recently to understand the experimen¬ 
tal observation of dipolar collisions of bosonic Feshbach 
Er 2 molecules in a ID optical lattice in a tilted magnetic 
field [3. 

The paper is structured as follows. In Section II, we 
develop the theoretical formalism for collision of parti¬ 
cles in a tilted field and confined geometry where the 
appropriate basis set for partial wave expansion is intro¬ 
duced. In Section III, we show how the collisional rate 
coefficients and their threshold behaviours are strongly 
affected by tilted fields revealing the complexity of the 
mechanism. The role of the confinement trap is also ex¬ 
plored and could be used to reduce the kinetic energy of 
the particles. Finally, we conclude in Section IV. 


II. THEORETICAL FORMALISM 

Ultracold dipolar collisions in a tilted field have 
been already studied in the past including microwave 
fields [36H38I ] but without confinement, in crossed elec¬ 
tric and magnetic fields or considering strong ID 

confinement such that the particles are bound to collide 
in pure 2D |4l|- By pure 2D we mean that the char¬ 
acteristic strength of the particles confinement is much 
stronger than the characteristic strength of their interac¬ 
tion (the dipole-dipole interaction here). However, this 
regime of pure 2D collisions is not yet reached in ongo¬ 
ing experiments as the required confinement strength is 
too strong. Instead, quasi-2D collisions occur. The par¬ 
ticles start to collide at large distances in pure 2D but 
there is a point as they approach each other where the 
increasing magnitude of their interaction gets much big¬ 
ger than the confinement strength. The particles do not 
feel anymore the presence of a 2D confinement and col¬ 
lide as if they were in a non-confined space. Therefore 
to reproduce the conditions of ongoing experiments, we 
describe the quasi-2D collisions of two ultracold parti¬ 
cles of mass mi, m 2 carrying tilted dipole moments and 
trapped in a ID optical lattice of arbitrary confinement 
strength. We assume that the particles cannot hop from 
one potential well to another so that we approximate 
a well by a harmonic oscillator for particle 1 and 2, 
Vho = 1/2 mi w 2 z\ + 1/2 m 2 w 2 z|. The angular fre¬ 
quency w = governs the strength of the confinement. 
The particles are initially in a given state of the harmonic 
oscillator 771,772 of energy e ni = hi/(ni + 1/2), e n2 = 
hv(ri2 + 1 / 2 ). 

The electric or magnetic fields make an angle 9e,b and 
<Pe,b with the z and x axis respectively as depicted in 
Fig. [2j-a. The classical dipole approximation has been 
shown convenient for modelling electric dipoles [34] and 
magnetic dipoles [19] interactions in ultracold dipolar 
gases. It is an appropriate way to avoid the inclusion 
of the particles internal structure in the collisional for¬ 
malism, thus sparing computational effort. No Stark nor 



FIG. 2: (Color online) a) Spherical angle coordinates 8e,b 
relative to the 2 axis and <£e,b relative to the x axis of the 
tilted fields E and B. In the study, the numerical results 
correspond to a field in the xOz plane with ifE,B = 0. b) 
Spherical coordinates (r, 9, ip) of the relative coordinate vector 
f describing the dipole-dipole collisional motion. 

Zeeman term appears in our formalism then. 

We distinguish three types of collisional processes: (i) 
Elastic processes for which the molecules remain in the 
same external state of the harmonic oscillator after the 
collision: (n'^n^) = (771,772); (ii): Inelastic processes 
for which the molecules change their external state of 
the harmonic oscillator (n'i,n' 2 ) ^ ( 771 , 772 ) (note that a 
change of the internal states is not possible as the inter¬ 
nal structure of the particles is not treated); (iii) Loss 
processes due for instance to chemical reactions occur¬ 
ring at ultralow energy, leading to products with high 
kinetic energy which are expelled from the trap. 


A. Quasi-2D collisions in parallel field 

We shall briefly recall the formalism for fields parallel 
to the quantization axis when 9e,b= 0 . It is presented in 
more details in Ref. [28J. First, it is convenient to trans¬ 
form the motion of the individual particles 1 and 2 with 
position coordinates r*i,r 2 of masses mi, m 2 in external 
states ?7 i, 772, into the motion of two effective particles 
described by the center-of-mass and relative coordinates 
R, r of masses rntot , 777 re d in external states TV, 77 . For 
confinements modelled by harmonic oscillators, it can be 
shown then that the center-of-mass coordinate R is de¬ 
coupled from the relative coordinate r. 

The potential energy term for the relative coordinate 
f is given by a van der Waals interaction 

V vdW = -C 6 /r 6 , ( 1 ) 

the confinement interaction given by a harmonic oscilla¬ 
tor for the relative motion along 2 

= 2 ^+ed a7 z , (2) 
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and a dipole-dipole interaction composed of an electric 
and magnetic term 


V M = V& + V& (3) 

where 

V dd m = —X^ m (l — 3 cos 2 9)/r 3 (4) 

with for electric dipoles 


X e = d/\/4n£ 0 (5) 

and for magnetic dipoles 

X m = p /\/ (47r//i 0 ). (6) 

The coordinate r is expressed in spherical coordinates 
r = (r, 9, <p) when V vd w and V dd m are dominant (Fig. [2} 
b), or in cylindrical coordinates r = ( p : z,tp ) when Vho is 
dominant at very large distances. 

We decompose the total wavefunction ijj(r) of the rel¬ 
ative motion into a basis set of spherical harmonics 
Y, m! ((9,vj) = (f\l,mi) 

V>(u #, ¥>) = ^ E ( r ) Y i mi V) (7) 

t,mi 


r, and we propagate the log-derivative of the wavefunc¬ 
tion [H, [ 43 . The boundary condition at short distance, 
where the propagation of the wavefunction is started, is 
set up by a tunable, diagonal log-derivative matrix given 
in Ref. [34] for which we can control the amount of loss 
in this short-range region. For the boundary condition at 
large distance in the asymptotic region where the prop¬ 
agation of the wavefunction is ended, we use the asymp¬ 
totic form of the cylindrical wavefunction which is a linear 
combination of regular and irregular Bessel functions of 
the cylindrical problem. Using a transformation matrix 
from cylindrical to spherical coordinates [28} , we obtain 
the asymptotic form of the wavefunction from which we 
deduce the K , S , and T matrices using the expression of 
the log-derivative matrix in spherical coordinates. The S 
matrix in the center-of-mass and relative coordinates is 
then expressed back into the individual coordinates [28} 
yielding the elastic, inelastic and loss cross sections and 
rate coefficients for two particles starting in a given initial 
state (ni, 712 ). 


B. Quasi-2D collisions in an arbitrary field 

In the spherical harmonics basis set 


with —l<mi< +1. In this basis set, the expression of As depicted in Fig. [2[ when the fields are tilted, the 
the potential energy terms become in bra-ket notation expression of the dipole-dipole interactions become 


(l mi\V vd w\l' m'i) 



X dmi -m'i 


( 8 ) 


(l mi\Vho\l' m'i) = 

1 2 2 
2 m red ^ r 

3 


mi ,m£ 

j(-i) m! V21TTV21' + 1 


and 


x 2 


l 2 l'\ 

0 0 0 ) 


l 2 l 1 
—mi 0 m\ 



(9) 


(l mi\V dd m \l' m'i) 


Vzox 


2 

e,m 


j.3 


s mi ,m[ (-i) mi V2iTiVWTT 

112Wi2 1'W l 2 V \ 
oooj^oooj^ -mi 0 m\ ) ' 


( 10 ) 


The first expression is diagonal in mi and l while the two 
last expressions are diagonal in mi but couple different 
values of l. The expression for V f lj(> arises from the fact 
that z 2 = r 2 cos 2 8 = r 2 [4 + l]/3. 

The time-independent Schrodinger equation is solved 
for a fixed total energy E to t = £m +£n 2 +E coil- Equations 
(0) to m lead to a set of coupled differential equations. 
To solve this system of equations, we use a diabatic-by- 
sector method [28| which generates a set of adiabatic en¬ 
ergy curves as a function of the inter-particle distance 


y d T = -Kmi 1 - 3cos 2 [0 - e E , B ])/r 3 . (11) 

This implies a more complicated expression in the spher¬ 
ical harmonic basis set 


(lmi\V d T\l'm?> = - 


Vwx 2 e , 


(-l) mi v / 2T+I V2F + 1 

112 

E E E 

Pl =-1 P 2 =—1 p= — 2 


1 1 2 
Pi P2 -P 


\JAir/SYf 1 ( 9e,b, <Pe,b) Y[ 2 (9e,b, 1 Pe,b) 


l 2 l' 

0 0 0 


l 


V 


-mi —p m\ 


( 12 ) 


For the special case Qe,b = 0, pi = P 2 = p = 0 and we 
recover the case mj —m; = 0. For the special case 9 e ,b = 
7t/2, pi = ±1, P 2 = ±1, p = 0, ±2 then we have m[— to; = 
0, ±2. For any other 9 e ,b we have m[—mi = 0, ±1, ±2, 
increasing the number of coupled equations and leaving 
no good quantum numbers in Eq. m- 


In the tesseral harmonics basis set 

By properly symmetrizing the basis set of spherical 
harmonics, we can still recover a good quantum number. 
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#E,B = 0, <PE,B = 0 

9e,b = 7t/2, ^Pe,b = 0 

9e,b 7 ^ 0, 7t/2, <pe,b = 0 

f>E,B ^ 0 

f = +l 

TOi = 0 
fhi — 1 
fhi = 2 

fhi = ... 

y 

y 

y 

y 

e = +i 

TOi = 
0 ,2,... 

TOi = 
1,3,... 

y 

y 

l—l 

+ 

II 

■Uy 

TO; = 

0,1, 2, 3,... 

y 

£ = ±i 

mi = 

0 ,1, 2,3,... 

y 

H-1 

TOi = 1 
fhi =2 

fhi = 3 

fhi = ... 

n 

n 

n 

n 

e = -i 

TOi = 

2.4.. .. 

TO; = 

1.3.. .. 

y 

y 

£ = -i 

TO; = 

1 ,2,3,4... 

y 


TABLE I: Quantum numbers £, to; needed ( y ) or not (n) for different angle 9e,b, <Pe,b of the fields. Note that for identical 
and indistinguishable (same internal state) bosons, l is even: l = 0, 2,4,... For identical and indistinguishable fermions, l is odd: 
Z = 1,3, 5,.... For distinguishable or non-identical particles, both parities of l should be included: l = 0, 2,4,... and l = 1,3, 5,.... 
For the special cases 9e,b = 0 and 7r/2 j_ if we start with molecules in the ground state of the harmonic oscillator, mi should be 
odd for fermions and even for bosons |28j|. 


We use the following symmetrized spherical harmonics in In this new basis set, the potential energy matrix ele- 


kct notation 




<5i g + z <5_i g 

V?A 

x {|i,-H)+e(-i) |mi| |i, Hi>} (is) 


where fhi = \mi\ and A = 1 + 6m,,o- The new quantum 
number £ takes the values £ = ±1 ,while fhi = 0,1, 2, 3... 
when £ = +1 and fhi = 1, 2,3... when £ = — 1. This new 
basis set are often called the tesseral harmonics [35] . We 
note them Y^m,,(, with 

o,c=+r = ^ {yr m + (-1 r if} 

oc P ; mi (cos0) cos fhi(f 

Yl,rn l= o, £ =+l = OS9) 

o,€=-i = -^{Y-^-i-ir'Y™ 1 } 

oc P™ ! (cos0) sinm;^. (14) 


ments are given by 

{l,fhi,£\V vdW \l'ifh'^i') = — j- x H' 6m,,m[ 6i,i>, (15) 


C 6 


(/, mi , £| Vh 0 \l , rrii, f; ) — _——. 6m,,m', 

VA A' 


x 2 


l 2 l' 
0 0 0 


l 2 V 

fhi 0 —to] 


+ 6i,i / + 6i,i> £' (— l) m ‘ dm,, o Hi 
and after reducing Eq. d 


(16) 




V30X 2 


r 3 VA 


(ci ) 2 


1 1 2 

1 1 -2 


l=(-ir^TT^FTl(J ^ J) 

g-i 2 ^e,b _|_^' e i2¥>E ’ 


i 2 v 

fhi —2 —to] 


(H' + 


+ 2 c 0 ci 

l 2 V 


1 1 2 

0 1 -1 


(H' 

i 


2 V 

-2 j 

g z <Pe,b e i<pE,B 


mi -i -to] J-^( -to. -i U -i to; 


2 l’ 

-f 


T 6m, ,m\ 


(co) 


2 , 1 1 2\ 2 / 1 12 
0 0 0 j 2( ' 1 ^ i -1 1 0 


l 2 V 

fhi 0 —to] 


l 2 V 

fhi 0 to] 


6m' 0 Hi 


(17) 
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c p,=0,1 



{0E,B, WE,B 


0 ), 


(18) 


in£,m;. Equation (El is diagonal in £ if ipE = 0 and/or 
(fB = 0, and it reduces in this case to 


namely cq = cos 0e,b and c\ = sin^.s/v^- Equa¬ 
tion El is diagonal in and Eq. El is diagonal 


J 


^=(-ift^nv2HT (* o o) 


(Cl) 2 

+ 2 cq ci 


1 1 2 

1 1 -2 

1 1 2 

0 1 -1 


V 

-I 


i : 

fhi —2 —fh l 

l 2 l 1 
fhi —1 —fh 




l 


l' 


(co) 5 


112 
0 0 0 


— 2( Cl ) 5 


1 1 2 

-1 1 0 


—mi —1 m[ 

l 2 V 

fhi 0 — mj 


+£'(-ir ; 


l 2 V 

fhi —2 fh'i 

l 2 l' 
fhi —1 mj 

2 / 


rn; 0 TO; 


• (19) 


We can see that we recover a good quantum number £ 
in this formalism which was not the case in the non- 
symmetrized spherical harmonics basis set. When <pe = 
<Pb 7 ^ 0 , or when a unique electric or magnetic field is 
used where lpe / 0 or / 0 , one can always recover the 
situation in Eq. El with good quantum numbers by a 
proper rotation of the x and y axis due to the cylindrical 
symmetry of the pancakes. Only in the more general 
case including both electric and magnetic fields where 
<pe 7 ^ one °f the dipole-dipole expression is no more 

diagonal in £ and Eq. El has to be used instead. To 
simplify the study in the following, we will consider the 
case where <pe,b = 0 when the dipoles are only tilted 
in the xOz plane. Furthermore we will consider the case 
where we only apply a unique field (electric in this study). 

In general, the quantum number fhi = 0 is always au¬ 
tomatically associated with the £ = +1 one as it is not 
defined for the £ = — 1 one. For the fhi > 0 values one has 
to use both £ = ±1 values. For the special case of paral¬ 
lel field 9e,b = 0, fhi is a good quantum number, C\ = 0 
and Eq. El reduces to Eq. El- Still for the parallel 
case for fhi 0, the £ = — 1 contribution is identical to 
the £ = +1 one so that the total contribution is twice the 
£ = +1 one. For fhi = 0, only the £ = +1 contribution 
is needed. Therefore, only the value of £ = +1 is needed 
in the parallel field case for all fhi. The contribution of 
the different quantum numbers needed for different tilted 
configurations are summarized in Tabic [T| 

The suitability of the tesseral harmonics can be qual¬ 
itatively understood from their spatial shape. If Lp = 0 
or 7 r which corresponds to the xOz plane (see Fig©, 
Yi t vanishes there for any values of fhi since it is 
proportional to sin(see Eq. El)- We can therefore 
associate the £ = — 1 quantum number with an out-of- 
plane motion, excluding the collision in the xOz plane. If 


<Pe,b = 0, the electric or magnetic field is applied in the 
xOz plane only. As the dipoles are exclusively pointing 
in the xOz plane, the collisions of dipoles in the y direc¬ 
tion start always side-by-side at long-range since there is 
no component of the dipole moment in the y direction. 
We thus expect that the manifold of the £ = — 1 curves 
is always more repulsive than its £ = +1 counterpart due 
to the side-by-side, repulsive dipolar approach. 

The tesseral harmonics representation of the partial 
waves is therefore a more appropriate basis set in the 
general case of collisions in arbitrary tilted electric and/or 
magnetic fields. Note that the present formalism is also 
adapted for collisions of particles in free space in crossed 
electric and magnetic fields, where one field (for example 
the magnetic field) is chosen as the quantization axis and 
the other (for example the electric field) is the tilted one. 
For example, the study done in Ref. j40| using spherical 
harmonics could be adapted using tesseral harmonics. 

III. RESULTS 

In the following, we study the collision of two indis¬ 
tinguishable, electric dipolar molecules of KRb in their 
absolute internal ground state for which the “classi¬ 
cal dipoles” assumption is a good description. We use 
40 K 8 'Rb for a fermionic system example and 41 K 8 'Rb 
for a bosonic one. We impose here ifiE = 0. They possess 
a permanent electric dipole moment of d p = 0.57 D Q. 
The electric field can therefore induce a dipole moment 
d up to d v 0. The results presented here will also be 
similar for collisions of magnetic dipolar particles with 
strong magnetic dipole moment, tuned by a tilted mag¬ 
netic field B 1 as performed in Ref. M- The confinement 
in the z direction is described by a harmonic oscillator of 
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FIG. 3: (Color online) Adiabatic energy curves as a function of the intermolecular distance r between two fermionic 40 K 87 Rb 
molecules (left panels) and two bosonic 41 K 87 Rb molecules (right panels) for a given induced dipole d = 0.2 D and tilt angle Be- 
The curves of same color are coupled among each other. Only the curves corresponding to £ = +1 are shown. For the £ = — 1 
curves, see Appendix A. Up and down arrows indicate the variation trend of the curves when d increases. Each adiabatic energy 
curve correlates to a state n of the harmonic oscillator indicated at the right of each panel. 


frequency v = 20 kHz which is a typical value employed 
in experiments. We assume that all the particles are in 
the ground state of the harmonic oscillator n\ = ri 2 = 0 
before the collision, which is equivalent to the relative 
n = 0 quantum number [28j . We study those collisions by 


varying different parameters such as the collision energy, 
the tilted angle and the confinement frequency. The loss 
processes are described by a full loss condition at short 
range given in Ref. [Tl]. The full loss condition corre¬ 
sponds to either a chemical reaction with full probability 
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at short range if the system is reactive, or if not to a 
possible “sticky” rate condition [45j where the two par¬ 
ticles stick together for a sufficient amount of time and a 
third one has the time to destroy the two-body complex 
equivalent to loss of particles. Although this mechanism 
has yet to be confirmed and observed in experiments, we 
include this possibility as well. We then consider the elas¬ 
tic, inelastic and loss rate coefficients for the four lowest 
values of to; = 0,1, 2,3 which are all mixed in a general 
tilted field 0<(?£<7r/2. In a parallel field 9e = 0, 
none of the to; are mixed and the corresponding rates 
for each fhi are summed altogether. In a perpendicular 
field 0e = 7t/2, the rates are calculated for the even value 
components fhi = 0, 2 and the odd value ones fhi = 1,3. 
For the bosonic case, the values of the mixed l are taken 
from l = 0 to l = 80 by steps of two, and for the fermionic 
case, from l = 1 to l = 79 by steps of two. The fact that 
we start with molecules in the ground state of the har¬ 
monic oscillator implies that for the special cases 9e = 0 
and 7r/2. to; should be odd for fermions and even for 
bosons 1281. 


A. Interactions and adiabatic energy curves 

The adiabatic energy curves are shown in Fig.[3]for the 
fermionic and the bosonic KRb molecules, for an induced 
dipole moment of d = 0.2 D. If there is a presence of a 
barrier relative to a given collision energy as r decreases 
and if this barrier increases when d increases (indicated 
by an arrow pointing upward in Fig. [3]), we say that the 
corresponding curve is protective against possible short- 
range loss. In the absence of a barrier or if the barrier 
decreases (indicated by an arrow pointing downward), we 
say that the corresponding curve is non-protective. The 
adiabatic energy curves are presented only for the quan¬ 
tum number £ = +1. For the £ = —1 quantum number, 
the curves are equal to their £ = +1 counterpart in the 
vanishing dipole moment limit, recovering the isotropic 
character of the long-range van der Waals interaction. 
For larger d, they are more repulsive due to the fact that 
the £ = — 1 manifold corresponds to side-by-side dipo¬ 
lar repulsive collisions as mentioned earlier. The £ = — 1 
curves are shown in Appendix A. 

The fermionic case for 9e = 0 is shown in the top 
left panel of Fig. [3] The thick black solid lines (resp. 
thin blue solid lines, thick red dashed lines, thin green 
dashed lines) represent the (unmixed) values of fhi = 1 
(resp. fhi = 3, to; = 0, fhi =2). The lowest curve of 
the n = 0 harmonic oscillator ground state correlates to 
the protective, repulsive side-by-side to; = 1 curve. For 
indistinguishable fermionic particles, we recall that the 
scattering takes place for odd l partial waves. The lowest 
curve l = 1 features a p- wave barrier. There is an actual 
crossing between the lowest curves with to; = 1 and fhi = 
0 since the fhi components do not couple to each other. 
The middle panel shows the case for 9e = 7t/ 4. The black 
solid lines represent the curves with mixed values to; = 


0,1,2, 3. Since all to; are coupled, the actual crossing 
at 9 e = 0 becomes an avoided crossing. And the lowest 
n = 0 curve correlates to a non-protective curve. Finally, 
the bottom panel shows the case for 9e = 7t/2. The 
black solid lines represent the curves with mixed values 
fhi = 1,3 and the red dashed lines the curves with mixed 
values fh; = 0, 2. Both series of curves do not couple to 
each other. The lowest n = 0 curve still correlates to 
the non-protective curve. So for particles in the lowest 
state, we can see how the adiabatic energy curve changes 
from a protective character coming from of a repulsive 
side-by-side interaction approach when 9e = 0 to a non- 
protective one coming from of an attractive head-to-tail 
interaction approach when the dipoles are tilted by 9e = 
tt/2. 

We get similar results in Fig. [3] for the bosonic sym¬ 
metry. For indistinguishable bosonic particles, the scat¬ 
tering takes place in even l partial waves. The lowest 
l = 0 curve is barrier less, in contrast with the fermionic 
particles. Now the lowest n = 0 curve correlates to the 
barrierless, non-protective curve at short range for all 9e 
even 9e = 0. As 9e increases from 0 to 7r/2, the fhi = 2 
curve pushes slightly downwards the fhi = 0 curve so 
that the lowest n = 0 curve get slightly more attractive. 

The behaviour of these adiabatic energy curves have 
direct consequences on the behaviour of the collisional 
rate coefficients. This is presented below. 


B. Collisions and rate coefficients 

1. Rate coefficients versus the collision energy. Threshold 

laws 

We show in Fig. [4]the rate coefficients for the fermionic 
and the bosonic case as a function of the collision energy 
for d = 0.2 D, starting with two particles in n = 0 rela¬ 
tive motional state. The top, middle and bottom panels 
corresponds to 9e = 0,7r/4,7r/2 respectively. The to¬ 
tal rates are reported as a thick solid line for each pro¬ 
cesses. One can see on Fig. [3] that when the collision 
energy is increased from the n = 0 threshold, different 
harmonic oscillator states become energetically open, re¬ 
sulting in sudden peaks in the inelastic rate coefficients 
(blue curves). In addition to the usual elastic rate coeffi¬ 
cient (red curves), we also have reactive rate coefficients 
(black curves) corresponding to the full loss condition at 
short range. 

For the fermionic particles for the 9e — 0 case, we 
see that reactive rates are suppressed compared to elas¬ 
tic one at ultralow energies. This has already been ex¬ 
plained theoretically and experimentally [26j, (28|, 123] ■ As 
the lowest n = 0 state curve correlates to a repulsive pro¬ 
tective barrier curve when r decreases, the probability of 
the two particles to be close to each other is low, pre¬ 
venting short-range loss to take place. The inelastic rate 
coefficient emerges at a collision energy of E c ~ 2 pK 
corresponding to the n = 2 threshold opening. The 
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FIG. 4: (Color online) Elastic (e), inelastic (i) and reactive (r) rate coefficients (thick lines) for fermionic (left panels) and 
bosonic (right panels) KRb + KRb collisions as a function of the collision energy at d=0.2 D and v = 20 kHz. Top panels: 
d E = 0, the thin solid lines are for fhi = 1 and £ = ±1 (to; = 0 and £ = +1), the thin dashed lines are for fhi = 3 and £ = ±1 
(fhi = 2 and £ = ±1) for fermions (bosons). Middle panels: 6e = 7r/4, the thin solid lines are for £ = +1, fhi = 0,1,2,3, 
the thin dashed lines are for £ = —1, fhi = 1,2,3. Bottom panels: 6e = tt/ 2, the thin solid lines are for £ = +1, fhi = 1,3 
(fhi = 0, 2), the thin dashed lines are for £ a* — 1, fhi = 1, 3 ( fhi = 2) for fermions (bosons). 


n = 1 threshold is not open here at E c ~ 1 /xK since 
the fhi = 0 and 2 (red and green curves) do not mix 
with the fhi = 1 and 3 (black and blue curves). When 
open above E c ~ 2 nK , the inelastic rates are bigger than 


the reactive ones and can easily amount or overcome the 
value of the elastic rates. For 0e = 7t/ 4, the reactive 
rates are bigger than the clastic ones at ultralow energies. 
This is explained again with the corresponding adiabatic 
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energy curve where the n = 0 state correlates to the at¬ 
tractive non-protective barrier due to the coupling of the 
fhi = 0 curve with the fhi = 1 one and the presence of 
the avoided crossing. The inelastic rate shows up now at 
E c ~ 1 /iK since the n = 1 threshold is now allowed and 
coupled with the n = 0 one. For 9e = 7t/ 2, we recover 
the head-to-tail collision with large reactive rate coeffi¬ 
cient. Due to the mixing of m; = 1,3 only, the n = 1 
threshold is not open here too and only the n = 2 opens 
up at E c ~ 2 /tK. 

For the bosonic particles, as the collision takes place 
on a barrier less l = 0 curve, we see that the magni¬ 
tude of the rates is bigger than the corresponding one for 
the fermions. For 6e = 0, elastic rates are comparable 
with the reactive ones (either smaller or bigger depend¬ 
ing on the collision energy) but there is no protection 
against collisions here. As 9 e increases for 9e = 7r/4 and 
7t/2, the reactive rate increases as the lowest adiabatic 
energy curve decreases and becomes the dominant rate 
coefficient. For the same reasons than those presented 
above for the fermions, the inelastic rates show up as the 
n = 2 threshold opens up at E c ~ 2 fi K for 9e = 0 and 
9e = 7t/2. For 9e = 7t/4, the n= 1 threshold also opens 
up at E c ~ 1 /tK. 

For all 9e ^ 0 plots, the contribution of the £ = — 1 
curves (thin dashed lines) is marginal for an induced 
dipole moment of d = 0.2 D. This is due to the large re¬ 
pulsive character of the corresponding adiabatic energy 
curves which correspond to side-by-side repulsive dipolar 
collision in the y direction, shown in Appendix A. 

For the inelastic collisions, there are striking differences 
at the threshold opening of the inelastic state depending 
on the field angle. This can be explained by the differ¬ 
ent behaviour of the quantum threshold laws studied in 
Ref. [46}. In this reference the threshold laws were de¬ 
rived for inelastic relaxation for (ki —> 0) -C kf where ki 
and kf represent the initial and final wavevector. Here 
we consider that the threshold laws for excitation pro¬ 
cesses are the same than the relaxation ones by replacing 
ki by kf and the initial mi i by the final mi f [47} . Then 
the inelastic relaxation laws given by Eq. (11) of Ref. [ib} 
where ki -C kf, for one or both of rrii i and mi f equal to 
zero, 


n 'i f 

An. relax. 


k. 


2 mi 


E c 


( 20 ) 


In 2 ki In 2 \J2m red E c 
should translate explicitly for the excitation processes to 

, 2mj , 


A. 


k " f 


(E c - E th ) m ‘f 


In kf In 2 y /2m ied (E c - E t h)/h 2 


■ ( 21 ) 


where kf -C ki, according to the behavior of the differ¬ 
ent elements of Eq. 10 of Ref. [4b[. E th corresponds to 
the excitation threshold energy. For both and mi f 
different than zero, the inelastic rate from Ref. [ibj] 


A 


in. relax. 


OC 


oc E c 


( 22 ) 


should translate to 


A 


mi,—>mi f 2mi f 

m. exc. ^ Kf 


oc(E c -E th ) mi f. 


(23) 


In our study, for the inelastic excitation rate n = 0 —> 
n = 2 at the E c ~ 2 /tK threshold for the specific cases 
9e — 0 and 9e — 7t/2, we found a threshold law of 
Annexe. K ( E c ~ -Eth) for the fermions and Annexe. a 
In -2 y / 2m re d(-E' c — E t if)/h 2 for the bosons, in agreement 
with Eq. ( 1251 ) and Eq. (Ell) , where E t h = e ra = 2 - At 
E c ~ 1 /iK, the inelastic excitation rate n = 0 —> n = 1 
for 9e = 7t/4 for the fermionic (resp. bosonic) case, gets 
the behaviour of the one for the bosonic (fermionic) case 
at 9e = 0 or 7t/2, with a sharper (rounder) shape. This is 
due to a change of parity in fhi in the inelastic transition. 
For the fermionic (bosonic) case, the inelastic transition 
n = 0 —> n = 1 corresponds to fhi = 1 —f fhi = 0 
(fhi = 0 —> fhi = 1 ). For the bosonic case for the 
n = 0 —> n = 1 inelastic transition, we found a threshold 
law of AdAexc. oc (E c - E th ) hr 2 ^2m red (E c - E th )/H 2 
(round shape) where E t h = e n =i now. For the 
fermionic case, we found a threshold law of Annexe. <x 
In -2 ^2m red (E c — E t h)/h 2 (sharp shape). Again this is 
in agreement with Eq. m- Finally, for the elastic and 
reactive collisions, the quantum threshold laws are found 
to be P e i oc E 2 and p le oc E c for the fermions in agree¬ 
ment with Eq. (12) and (14) determined in Ref. [46( for 
dipolar collisions in quasi-2D. For the bosons we find 
Pei oc p re oc ln~ 2 ^2m red E c /h 2 in agreement with Eq. (9) 
and (13) of the same reference. 


2. Sensitivity of the rate coefficients versus the tilted field 

angle 

The effect of the tilted field on the collision is directly 
seen in Fig. [5] for fermions (top panel) and bosons (bot¬ 
tom panel). This is plotted at a fixed collision energy of 
E c = 1 /iK and induced dipole of d = 0.2 D. 

There is no fermionic and bosonic inelastic rate at 9e = 
0 and 9e = 7t/ 2 at E c = 1/rK. Couplings between fhi = 1, 
fhi = 3 and fhi = 0, fhi — 2 curves are not allowed at 
these two angles so that the transition n = 0 — >-n=lis 
forbidden there. Then the inelastic rates rise from 9e = 0 
to 9e = 7t/ 4 due to the turning on of the couplings. At 
9e = 7t/ 4 there is a maximal coupling between the fhi 
components: (i) a maximal coupling for A fhi = 2 due to 
the presence of the sin 2 9e = 1/2 term in Eq. (THU) : (ii) 
a maximal coupling for Am; = 1 due to the presence of 
the cos 9e sin 9e = 1/2 term. Thus the inelastic rates 
reach their maximal value at 9e = 7t/ 4. Conversely, the 
couplings turn off as 9e passes from 7 t/ 4 to 7r/2, and then 
the inelastic rates shut off. 

The fermionic reactive rate increases as we pass from 
a side-by-side approach 9e = 0 to a head-to-tail one 
9e = 7r/2. The continuous transition seen in the rate as 
9e increases can be understood qualitatively by the con¬ 
tinuous increase of the A fhi = 1 coupling. At 9e = 0 the 
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lowest curve of symmetry fhi = 1 connects to the protec¬ 
tive barrier and leads to a rate Po, while at Be = 7t/2 this 
curve connects to the non-protective barrier that leads to 
a rate P v / 2 » Po- Now, at Be = 0 the lowest curve of 
symmetry fhi = 0 connects to the non-protective barrier, 
the same one that leads to the rate p^m. Therefore as 
the coupling between the two symmetries increases from 
Be = 0 to 7 t/ 4, the reactive rate is a combination of 
Po and 13^ / 2 with coefficients that decrease the contribu¬ 
tion of Po and increase the one of P„/ 2 , hence increasing 
the total reactive rate. From Be = 7r/2 to 7r/4 the re¬ 
verse argument holds since now at Be = 7r/ 2 the lowest 
curve of symmetry to; = 1 (to; = 0) connects to the 
non-protective (protective) barrier. The reactive rate is 
again a combination of 0^/2 and po- But as the coupling 
increases from Be = 7t/ 2 to 7r/4, the coefficient of P w / 2 
decreases while the one of Po increases, hence decreas¬ 
ing the total reactive rates and connecting to the trend 
between 0 E = 0 and 7r/4. 

The sensitivity of the fermionic reactive rate with the 
field angle for the KRb system is quite strong, since a 
small change of Be = 7r/10 (= 18°) gives rise to an order 
of magnitude increase in the rate. Therefore in exper¬ 
iments of fermionic KRb molecules, it is important for 
the electric field to be quite parallel to the ID optical 
lattice confinement axis to avoid additional losses due to 
slight tilted angles. Note however that the range of this 
strong sensitivity of the rates may vary from one system 
to another. 

In contrast for bosons, there is a very slight dependence 
of the reactive rate with Be- This is due to the slight 
downward pushing of the to; = 2 curve to the fhi = 0 
ones when Be increases as can be barely seen on Fig. [3] 
The reactive rate value is anyway larger than the elastic 
and inelastic one. 

Additionally, it is interesting to analyze the effect of 
the £ component. For fermions, the £ = — 1 rate coef¬ 
ficient turns out to be the same than the £ = +1 one 
at Be = 0 because Eq. (1191) reduces to Eq. (1101) which 
is an expression independent of the £ number. When 
departing from the angle Be = 0, the £ = +1 rate coeffi¬ 
cients are dominant compared to the £ = — 1 ones. As a 
good approximation, one can neglect the contribution of 
the latter for large values of Be (note that as d increases, 
this is less true though, see comment in Appendix A). For 
bosons at Be = 0, the £ = +1 term contains the to; = 0 
component responsible for a barierless collision while in¬ 
trinsically the £ = — 1 one does not. Therefore the latter 
case does not yield a higher rate than the £ = +1 case, 
then it is a good approximation to neglect the £ = — 1 
component for bosons for all angles Be- 


3. Rate coefficients versus the confinement 

The effect of the confinement strength is shown in 
Fig. [G] for a fixed collision energy of E c = 500 nK, in¬ 
duced dipole of d = 0.2 D and field angle Be = 7r/10 



0 E (71) 



FIG. 5: (Color online) Elastic (e), inelastic (i) and reactive 
(r) rate coefficients (thick lines) for fermionic (top panel) and 
bosonic (bottom panel) KRb + KRb collisions as a function 
of the tilt angle Be, at E c =1 fiK, d = 0.2 D and v = 20 kHz. 
The partial rates are given by thin solid lines (£ = +1, fhi = 
0,1, 2, 3) and thin dashed lines (£ = —1, m; = 1, 2, 3). 


(18°). At v = 20 kHz, the first threshold m = 0,772 = 1 
is located at an energy of ~ 1/rK above the energy of 
the initial state (ni = 0,77-2 = 0). Therefore a collision 
energy of 500 nK is not sufficient to open up inelastic 
collisions and there is no inelastic rate. When the con¬ 
finement decreases, there is a given v for which the first 
inelastic transition becomes open, when the first thresh¬ 
old (?7i = 0,772 = 1) energy amounts the value of the 
collision energy. This is the case here at v = 10.4 kHz. 
These different values of v are indicated with an arrow 
along with the different threshold openings ( 771 , 772 ). At 
v = 10.4 kHz, the first inelastic threshold opening (first 
arrow from the right) for the fermions is quite strong, 
recalling the sharp one seen in Fig. [4] for Be = 7r/4. 
The first one for the bosons is rather weak and smooth 
as also seen in Fig. [4j The second opening (second ar¬ 
row from the right) is now smooth for the fermions and 
sharp for the bosons. And so forth, the successive open¬ 
ings alternate between sharp and smooth patterns. As 
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v (kHz) 


FIG. 6: (Color online) Elastic, inelastic and reactive rate coef¬ 
ficients for fermionic (top panel) and bosonic (bottom panel) 
KRb + KRb collisions as a function of the tilt angle 9e, at 
E c =500 nK, d=0.2 D and v = 20 kHz. The numbers in the 
brackets represent the harmonic states (ni,ri 2 ) of particles 1 
and 2 and the arrows the corresponding energy thresholds. 
The inset in the top panel sketches the excitation inelastic 
process (0,0) —» (0,1) (see text for details). 


explained earlier, the smooth openings correspond to a 
fhi = 0 — i fhi = odd transition while the sharp openings 
correspond to a fhi = 1 —> fhi = even transition. 

An interesting feature is seen for the fermionic sys¬ 
tem. At v = 9 kHz for example, the inelastic rate 
(ni = 0 ,ri 2 = 0) —> (m = 0 ,ri 2 = 1) reaches the 
value of the elastic rate (~ 2 10 -6 cm 2 s -1 ) while the 
reactive rate is about 3.3 times smaller. The inelastic 
process here is an excitation inelastic process (see the 
inset of Fig. [6]). The total energy is set before the col¬ 
lision by E tot = £ 0 + e 0 + E c = 0.932/iK. We also have 
£o = 0.216/iK and £\ = 0.648/iK. After the collision, the 
total energy is given by E tot = e 0 + £\ + e[ where El 
is the finall kinetic energy. The conservation of the total 
energy implies = E c — (e-^ — e 0 ) = E c — hv = 500 
nK—432 nK= 68 nK. This means that for an excitation 
inelastic process, the particles after the collision have a 


smaller kinetic energy than before the collision and there¬ 
fore slowed down by this mechanism. Loss of molecules 
can happen but according to the rates, for each loss of one 
pair, three pairs get slowed down from 500 nK to 68 nK. 
Note that the final kinetic energy depends on the choice 
of the frequency trap (see inset): its value is even smaller 
when the frequency is closer to the frequency that shuts 
off the inelastic transition, here v = 10.4 kHz. Therefore, 
particles can be slowed down to arbitrary small kinetic 
values, as far as the reactive rates remain smaller than 
the inelastic ones. Of course the reverse relaxation in¬ 
elastic process (m = 0 ,n 2 = 1) — > {n\ = 0 ,ri 2 = 0) 
can also happen after the particles have been excited, 
restoring back the 500 nK kinetic energy to the parti¬ 
cles. But this can be prevented by tilting the electric 
field back to the parallel case 6e = 0. In such cases, the 
inelastic transition is forbidden as discussed previously, 
leaving the particles in the harmonic states 0 and 1 with 
small kinetic energy. Another possibility is to find a way 
to remove directly the particles in the harmonic state 1 
leaving only the one in the ground state 0 with small 
kinetic energy. 


IV. CONCLUSION 


By implementing tesseral harmonics in place of spher¬ 
ical harmonics in the collisional formalism of two ultra¬ 
cold tilted dipolar particles in confined space, we showed 
that we can recover a good quantum number £. This 
separates the overall problem into two sub-problems of 
smaller size even when a field is tilted. Inelastic and reac¬ 
tive rates show dramatical changes in a tilted field. This 
is due to additional couplings between fhi components 
when a tilt is applied. We also showed that fermionic 
dipolar particles can lose kinetic energy and slow down 
due to favourable trap excitation inelastic collision under 
an appropriate confinement strength of the ID lattice in 
a tilted field. Future works will investigate whether this 
mechanism can be efficient to even further cool down the 
particles in such configuration by taking into account the 
initial kinetic energy distribution of the particles for a 
given temperature and their rethermalization due to the 
elastic colliisons during this process. 
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Appendix A: £ = ±1 adiabatic energy curves 

Fig- Q presents the adiabatic energy curves for d = 
0.2 D and 9e = 7r/4 for both the fermionic and bosonic 
system and both quantum numbers £ = ±1. One can 
see that the £ = — 1 curves always corresponds to more 
repulsive curves than £ = +1 ones. This is due to the fact 
that fhi > 1 for £ = — 1 and it correlates to a l > 1 curve 
with a large centrifugal barrier. At weak d the £ = — 1 
manifold does not play an important role in the collision. 
Note that at large d , the l > 1 curve can turn attractive 
again and then the £ = —1 manifold can start to play a 
role in the collision |28j . 



FIG. 7: (Color online) Adiabatic energy curves as a function 
of the intermolecular distance r for d = 0.2 D, 0e = rr/4 
for the £ = +1 manifold (black solid lines) and the £ = — 1 
manifold (red dashed lines). Fermions: top panel, bosons: 
bottom panel. Each adiabatic energy curve correlates to a 
state n of the harmonic oscillator. 


[1] J. Doyle, B. Friedrich, R. Krems, and F. Masnou-Seeuws, [2] M. Baranov, Phys. Rep. 464, 71 (2008), ISSN 0370-1573. 
Eur. Phys. J. D 31, 149 (2004). [3] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and 





















































13 


T. Pfau, Rep. Prog. Phys. 72, 126401 (2009). 

[4] G. Quemener and P. S. Julienne, Chem. Rev. 112, 4949 
( 2012 ). 

[5] M. A. Baranov, M. Dalmonte, G. Pupillo, and P. Zoller, 
Chem. Rev. 112, 5012 (2012). 

[6] S. Kotochigova, Rep. Prog. Phys. 77, 093901 (2014). 

[7] K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe’er, 
B. Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, 

D. S. Jin, and J. Ye, Science 322, 231 (2008). 

[8] K. Aikawa, D. Akamatsu, M. Hayashi, K. Oasa, 
J. Kobayashi, P. Naidon, T. Kishimoto, M. Ueda, and 

S. Inouye, Phys. Rev. Lett. 105, 203001 (2010). 

[9] T. Takekoshi, L. Reichsollner, A. Schindewolf, J. M. Hut¬ 
son, C. R. Le Sueur, O. Dulieu, F. Ferlaino, R. Grimm, 
and H.-C. Nagerl, Phys. Rev. Lett. 113, 205301 (2014). 

[10] P. K. Molony, P. D. Gregory, Z. Ji, B. Lu, M. P. 
Koppinger, C. R. Le Sueur, C. L. Blackley, J. M. Hutson, 
and S. L. Cornish, Phys. Rev. Lett. 113, 255301 (2014). 

[11] J. W. Park, S. A. Will, and M. W. Zwierlein, Phys. Rev. 
Lett. 114, 205302 (2015). 

[12] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and 

T. Pfau, Phys. Rev. Lett. 94, 160401 (2005). 

[13] Q. Beaufils, R. Chicireanu, T. Zanon, B. Laburthe-Tolra, 

E. Marechal, L. Vernac, J.-C. Keller, and O. Gorceix, 
Phys. Rev. A 77, 061601 (2008). 

[14] B. Naylor, A. Reigue, E. Marechal, O. Gorceix, 
B. Laburthe-Tolra, and L. Vernac, Phys. Rev. A 91, 
011603 (2015). 

[15] M. Lu, N. Q. Burdick, S. H. Youn, and B. L. Lev, Phys. 
Rev. Lett. 107, 190401 (2011). 

[16] M. Lu, N. Q. Burdick, and B. L. Lev, Phys. Rev. Lett. 
108, 215301 (2012). 

[17] K. Aikawa, A. Frisch, M. Mark, S. Baier, A. Rietzler, 
R. Grimm, and F. Ferlaino, Phys. Rev. Lett. 108, 210401 
( 2012 ). 

[18] K. Aikawa, A. Frisch, M. Mark, S. Baier, R. Grimm, and 

F. Ferlaino, Phys. Rev. Lett. 112, 010404 (2014). 

[19] A. Frisch, M. Mark, K. Aikawa, S. Baier, R. Grimm, 
A. Petrov, S. Kotochigova, G. Quemener, M. Lepers, 
O. Dulieu, and F. Ferlaino, ArXiv 1504.04578. 

[20] B. K. Stuhl, M. T. Hummon, M. Yeo, G. Quemener, J. L. 
Bohn, and J. Ye, Nature 492, 396 (2012). 

[21] J. F. Barry, D. J. McCarron, E. B. Norrgard, M. H. Stei- 
necker, and D. Demille, Nature 512, 286 (2014). 

[22] M. Yeo, M. T. Hummon, A. L. Collopy, B. Yan, B. Hem- 
merling, E. Chae, J. M. Doyle, and J. Ye, Phys. Rev. 
Lett. 114, 223003 (2015). 

[23] A. L. Collopy, M. T. Hummon, M. Yeo, B. Yan, and 
J. Ye, New J. Phys. 17, 055008 (2015). 

[24] B. Pasquiou, A. Bayerle, S. M. Tzanova, S. Stellmer, 


J. Szczepkowski, M. Parigger, R. Grimm, and F. Schreck, 
Phys. Rev. A 88, 023601 (2013). 

[25] C. Ticknor, Phys. Rev. A 81. 042708 (2010). 

[26] G. Quemener and J. L. Bohn, Phys. Rev. A 81, 060701 

( 2010 ). 

[27] A. Micheli, Z. Idziaszek, G. Pupillo, M. A. Baranov, 
P. Zoller, and P. S. Julienne, Phys. Rev. Lett. 105, 
073202 (2010). 

[28] G. Quemener and J. L. Bohn, Phys. Rev. A 83, 012705 

( 2011 ). 

[29] M. H. G. de Miranda, A. Chotia, B. Neyenhuis, D. Wang, 

G. Quemener, S. Ospelkaus, J. Bohn, J. L. Ye, and D. S. 
Jin, Nature Physics 7, 502 (2011). 

[30] A. Chotia, B. Neyenhuis, S. A. Moses, B. Yan, J. P. 
Covey, M. Foss-Feig, A. M. Rey, D. S. Jin, and J. Ye, 
Phys. Rev. Lett. 108, 080405 (2012). 

[31] A. Simoni, S. Srinivasan, J.-M. Launay, K. Jachymski, 
Z. Idziaszek, and P. S. Julienne, New J. Phys. 17, 013020 
(2015). 

[32] B. Yan, S. A. Moses, B. Gadway, J. P. Covey, K. R. A. 
Hazzard, A. M. Rey, D. S. Jin, and J. Ye, Nature 501, 
521 (2013). 

[33] L. D. Carr, D. DeMille, R. V. Krems, and J. Ye, New J. 
Phys. 11, 055049 (2009). 

[34] G. Wang and G. Quemener, New J. Phys. 17, 035015 
(2015). 

[35] E. T. Whittaker and G. N. Watson, Cambridge Univer¬ 
sity Press (1990). 

[36] A. V. Avdeenkov, New J. Phys. 11, 055016 (2009). 

[37] A. V. Avdeenkov, Phys. Rev. A 86, 022707 (2012). 

[38] A. V. Avdeenkov, New J. Phys. 17, 045025 (2015). 

[39] E. Abrahamsson, T. V. Tscherbul, and R. V. Krems, J. 
Chem. Phys. 127, 044302 (2007). 

[40] G. Quemener and J. L. Bohn, Phys. Rev. A 88, 012706 
(2013). 

[41] C. Ticknor, Phys. Rev. A 84, 032702 (2011). 

[42] B. R. Johnson, J. Comp. Phys. 13, 445 (1973). 

[43] D. E. Manolopoulos, J. Chem. Phys. 85, 6425 (1986). 

[44] K.-K. Ni, S. Ospelkaus, D. Wang, G. Quemener, 
B. Neyenhuis, M. H. G. de Miranda, J. L. Bohn, D. S. 
Jin, and J. Ye, Nature 464, 1324 (2010). 

[45] M. Mayle, G. Quemener, B. P. Ruzic, and J. L. Bohn, 
Phys. Rev. A 87, 012709 (2013). 

[46] Z. Li, S. V. Alyabyshev, and R. V. Krems, Phys. Rev. 
Lett. 100, 073202 (2008). 

[47] H. R. Sadeghpour, J. L. Bohn, M. J. Cavagnero, B. D. 
Esry, I. I. Fabrikant, J. H. Macek, and A. R. P. Rau, J. 
Phys. B: At. Mol. Opt. Phys. 33, 93 (2000). 


